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Abstract. We explore a plausible mechanism that the hemispherical power asymmetry in 
the CMB is produced by the spatial variation of the primordial sound speed parameter. We 
suggest that in a generalized approach of the 6N formalism the local e-folding number may 
depend on some other primordial parameters besides the initial values of inflaton. Here the 6N 
formalism is extended by considering the effects of a spatially varying sound speed parameter 
caused by a super-Hubble perturbation of a light field. Using this generalized 6N formalism, 
we systematically calculate the asymmetric primordial spectrum in the model of multi-speed 
inflation by taking into account the constraints of primordial non-Gaussianities. We further 
discuss specific model constraints, and the corresponding asymmetry amplitudes are found to 
be scale-dependent, which can accommodate current observations of the power asymmetry at 
different length scales. 
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1 Introduction 

The hemispherical power asymmetry in the cosmic microwave background (CMB) was indi¬ 
cated by the WMAP experiment years ago [1-4]. Later the Planck collaboration also reported 
this feature in the CMB fluctuations [5, 6] and provided an independent measurement of this 
anomaly. Phenomenologically, such a power asymmetry can be modeled as a dipolar modu¬ 
lation of a statistically isotropic CMB sky in terms of temperature fluctuations in a specific 
direction of the unit vector p. Translated into the expression of the primordial power spec¬ 
trum, the modulation required to explain this asymmetry can be written as a spatially-varying 
power spectrum, 


P^{k,f) = P(;{k)[l+ 2A{k) p-{K-:xio)/xis] , (1.1) 

where P({k) is the isotropic spectrum, A{k) characterizes the amplitude of dipolar asymmetry, 
xis is the comoving distance to the last scattering surface and xq represents our current 
location. 

The amplitude of the dipolar asymmetry is given by A = 0.072 ± 0.022 for the CMB 
power with < 60 [5]. However, the asymmetry does not necessarily exist at smaller length 
scales. Particularly, the constraint from the Sloan Digital Sky Survey quasar sample [7] 
requires A < 0.0153 (99% C.L.) for the power asymmetry oriented in the direction of the 
CMB dipole in which the typical wavenumber is A: ~ lMpc“^. And analyzing CMB data only 
also indicates a small asymmetry amplitude at large i [8, 9]. Thus, any model that accounts 
for the CMB power asymmetry has to produce a strong scale dependence so that it can be in 
agreement with the above observational constraints. In order to explain the observed scale- 
dependent CMB anisotropy various theoretical models have been proposed in the literature 
[10-41]. 

As pointed out in [14], a single-field slow-roll inflation model cannot generate such an 
asymmetry without violating the constraints of homogeneity in the universe. The same pa¬ 
per also proposed a so-called Erickcek-Kamionkowski-Carroll (EKC) mechanism based on a 
curvaton model [42-45] and thus can explain this anomaly without violating the homogeneity 
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constraint. However, the original model is inconsistent with the quasar bound since the sig¬ 
nature is scale-independent. Also, the model leads to a large value of non-Gaussianity which 
has been ruled out by Planck [46]. Instead, a improved curvaton model in which the curvaton 
decay takes place after dark matter freezes out was studied [47]. 

When dealing with cosmological perturbations in a system involving multiple scalar 
fields, it turns out that the 6N formalism [48-52] is very powerful, particularly when describing 
nonlinear terms. The application of the 6N formalism in the study of the CMB power 
asymmetry has been carried out in the literature [14, 20, 25, 28, 30, 31, 34-36]. It was found 
in these studies that there may exist certain discrepancy if the nonlinear contributions to the 
quadrupole were not treated in a full analysis, as was pointed out in [36]. 

In spite of all those models, we suggested that the hemispherical power asymmetry may 
be generated by a modulation of primordial speed of sound parameter during inflation as 
studied in Ref. [29]. In the inflation model of K-essence field [53], the propagation of the 
field fluctuation was found to be governed by a speed of sound parameter c^. This parameter 
usually equals to unity when the kinetic term is of standard form. However, when the kinetic 
term is nontrivial, the sound speed parameter could be smaller than unity and yields a different 
squeezing process for the field fluctuation and thus modifies the primordial power spectrum 
of the curvature perturbation. Further, the model of multi-speed inflation is constructed by 
a collection of scalar fields with such nonstandard kinetic terms [54, 55]. Then, in this model 
there is one corresponding sound speed for each field, and therefore the fluctuations of these 
fields do not propagate synchronously. 

In the present work we extend the 5N calculation by considering the spatially varying 
primordial sound speed parameter of the inflaton field. Such a modulation of the sound speed 
can be obtained by the super-Hubble perturbation of a second light field, which is dubbed 
the EKC effect. This mechanism is implemented in the model of the multi-speed inflation. 
In Section H, we give a brief review of this model and the 5N formalism. In Section HI we 
illustrate the feasibility of our mechanism by giving a detailed calculation of the asymmetry 
amplitude A{k). Afterwards, in Section IV we consider a specific model and numerically 
examine the effect of the modulation of sound speed parameter. Finally a scale-dependent 
asymmetry amplitude is obtained to confront the observational data. Section V includes the 
conclusion with a discussion. 


2 A brief review of multi-speed inflation and 6N formalism 


Our starting point is based on the multi-speed inflation model, which was proposed in Refs. [54, 
55]. In a simple version of the model we consider that the universe is driven by two scalar 
fields minimally coupled to Einstein gravity at early times, with one being a canonical field 
X and the second being a K-essence field (j) which has a non-canonical kinetic term P{X,(j)) 
with X = — This K-essence type of the scalar field may arise from a low 

energy effective field description of string theory or other more fundamental theories [53]. For 
instance, the non-canonical kinetic term could be of the Dirac-Born-Infeld-like (DBI) form, 
and accordingly, the Lagrangian density can be written as 

C = P{X,cP,x)-ld^xd>^X-Vic^,x) , ( 2 . 1 ) 


with 


P{X,(^,X) 


1 


\/l+ /(</', 


( 2 . 2 ) 


2 






The coefficient f((j), x) inside the square root function is a warping factor for the scalar held 


(j). 

Here we assume that cf) plays the role of inflaton while y makes no contribution to 
the background dynamics^. In this model a distinct property is that the speed of sound 
parameters of those two helds can behave differently. For x, the sound speed is simply unity 
since its kinetic term is canonical. For (p, however, the square of the sound speed is given by 


2 ^ P,x 
P,x + 2XPxx 


l-2f{cP,x)X , 


(2.3) 


where Px denotes the partial derivative of P with respect to X. Note that Cg is time- 
dependent and affected by the variation of scalar held x- this model the sound speed 
parameter is expected to be related to the primordial non-Gaussianities via the following 
expression [59] 

Then one can observe that a small value of the sound speed parameter can lead to a large 
amount of primordial non-Gaussianity, while c* ~ 1 corresponds to /nl ^ 1- Hence, the 
model can be well constrained by cosmological observations such as the Planck data [60, 61]. 


2.1 Background dynamics 

To study the background dynamics, we vary the Lagrangian with respect to the metric and 
then obtain the energy-momentum tensor as follows, 


^ _ y) + PxV^0V^(/> . (2.5) 

Moreover, varying the Lagrangian with respect to the inhaton yields the following generalized 
Klein-Gordon equations. 


V^(PxV^(/)) + P^-P,^ = 0 , (2.6) 

where P,^ denotes the partial derivative of P with respect to the scalar p. 

Gonsidering a spatially hat Friedmann-Robertson-Walker (FRW) spacetime with its met¬ 
ric 

ds^ = —dt^ + a^{t)d-^ , (2.7) 

we can read the energy density and pressure of the above cosmological system from the 
energy-momentum stress tensor 


p = —-P + R, p = P-V. 2.8 

The equations of motion for the scalar helds now can be expressed as 

■ c ■ dP dV 

p + 3H(j) - -(j) - Cg—+ Cg—= 0 , (2.9) 

Cg op op 

^For instance, this model was applied in Refs. [56, 57] in order to realize a viable curvaton mechanism and 
was investigated in [58] to produce negative-valued local non-Gaussianities. 
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where the Hubble expanding rate H = aja has been introduced. 
Equation is given by 


H = - 


2^2 


{P + P) = - 


2M?,c, 


V- 


The second Friedman 
( 2 . 10 ) 


where Mpi is the reduced Planck mass, which is defined as l/\/8vrG. Moreover, we can 
introduce the following three slow roll parameters 


_ H _ e _ Cs 

’ ^ = jre' 


( 2 . 11 ) 


and they respectively characterize the variation of H, e and Cg in one Hubble time. To ensure 
a sufficiently long and stable inflationary phase, these parameters are demanded to be much 
smaller than 1 during inflation. In addition, combining Eqs. (2.10) and (2.11), we get the 
following relation 

^ = 2M2ec. ’ 

which is very useful in the analysis of the generalized 5N formalism in next section. 


2.2 5N formalism and the isotropic spectrum 

Having the background dynamics of multi-speed inflation, we review the basic idea of the 6N 
formalism in this subsection. Afterwards, we use this approach to calculate the power spec¬ 
trum of primordial curvature perturbation P(^{k). During inflation the regions of super-Hubble 
scales may be regarded as many “separate universes”. Each of these regions is approximately 
homogeneous and isotropic, and thus, the corresponding evolution can be locally described by 
an FRW one. However, due to quantum fluctuations, the primordial parameters, such as the 
initial value of inflaton, are different in these regions, which can lead to different expansion 
behavior. We may dehne N{x, t) as the local e-folding number from a spatially flat slice 
to a uniform-density slice at time t. Then the curvature perturbation of this region can be 
interpreted as the difference between N{x, t) and the background e-folding number, i.e. 


C (x, t) = N{x,t) — N(t) = 6N{x, t) , 


(2.13) 


which is known as the 5N formalism and is formulated as follows. 


6N 


dN 1 d^N 
dcpi * 2 d(j)id(j)j 


64>i6(t)j -|-... , 


(2.14) 


where a subscript j denotes the evaluation along the i-th unperturbed field trajectory. 

Here 5(pi corresponds to the i-th field fluctuations whose wave lengths are smaller than 
the current observable universe, i.e. k > Since in our model the background evolution 
is independent of Xj there is ^ = 0. In this case, the relation between N and (j) is simply 
determined by the background dynamics of inflation 


dN _ H 
d(j) ^ 


(2.15) 


To leading order in the expansion of Eq. (2.14), the power spectrum of curvature perturbation 
turns to be 


Pdk) 


k^ 

27r^ 


ON 

dtp 




(2.16) 
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As is well known, = ■^\54>k\'^ 
(2.15) and (2.12), we derive 


in a quasi-de Sitter spacetime. 


dN 


d(f) 


2My, 


Then using Eqs. 


(2.17) 


and thus, the isotropic power spectrum in our model is expressed as. 


^ ~ 87r2M2ec. 


(2.18) 


3 Asymmetric spectrum caused by the spatial varying sound speed 

In this section we firstly extend the 5N formalism by taking into account the effects of a 
varying speed of sound parameter. Then we use this generalized approach to analyze the 
generation of power asymmetry caused by the spatially varying sound speed. 


3.1 A generalized 5N formalism 

As was reviewed in the previous section, the regular 5N formalism considers only the direct 
dependence of the e-folding number N on the inflaton cj). However, there exist some primor¬ 
dial parameters that may vary with the background evolution during inflation. In our case, 
from Eqs. (2.9) and (2.10), one can observe that the primordial sound speed parameter Cg is 
such a parameter. Thus we suggest that besides the initial value of inflaton <^(x), Cg also influ¬ 
ences the local e-folding number through a generalized formalism: N{x, t) = N{(f)(x),Cg{x),t). 
Accordingly, the 6N formula in Eq. (2.14) can be extended as follows. 


, ON ^ 

^^ ^ ^ 2W * ^ dc,d4, 


(3.) 


up to second order. As before, the derivatives of N with respect to Cg are evaluated at the 
value of the background sound speed Cg. In the following, we use the extended 6N formalism 
to analyze the effects of a spatially varying sound speed parameter and apply this approach 
to calculate the corresponding asymmetric primordial power spectrum. 


3.2 The spatial variation of the sound speed 

Let us consider the situation that the sound speed Cg is varying spatially as follows. 


Cg (x, t) = Cs (t) Acg (x, t) . 


(3.2) 


According to Eq. (2.3), this modulation of the sound speed can be achieved by considering 
the long wavelength field perturbation Ax(x,t) through 


Acg{x,t) 


dcg 

dx 


Ax(x,t) 


(3.3) 


From Eqs. (2.3) and (2.4), we have 


dcg 

dx 



- 5 - 


c, 


(3.4) 










For Jnl < 1 (i-e. Cg ~ 
[29], which is given by 


1), the leading order of Eq. (3.3) reduces to the one developed in Ref. 


Cfi — Cc 


1 + 

35 / 


Different from the sub-Hubble perturbation 6x, here Ax(x, t) is a super-Hubble mode with 
specific direction and wavelength. According to the Grishchuk-Zel’dovich (GZ) effect [62], the 
very large scale mode Ax could modify the spectrum in the observable universe. 

Following the treatment in Ref. [36], we assume that Ax only exists at single wave 
number and treat it as a non-stochastic quantity. Thus it can be expressed as 


Ax(x, t) = B{t) cos(kL • x) , 

where ki, We further assume that we are located at xq, then there is 

ki • X = ki • xo ki ■ (x - xo) . 

Gonsidering the limit of k^, • (x — xq) <C 1, one can expand Eq. (3.5) to be 

Ax(x, t) = B{t)[ cos 6 — sin0 k/, • (x — xq)] , 


(3.5) 


(3.6) 


(3.7) 


with 0 = ki • Xq. As a result, the spatially varying part of the sound speed parameter takes 
the form of 

- -- 

ACs(x, t) = Acs(f) -|- ACs(x, t) = -;:^R(t) [cos 0 — sin0 k/, • (x — Xq)] . (3.8) 

ox 

3.3 The spatially varying spectrum 

Having the above consideration, now we turn to calculate the asymmetric spectrum. In the 
Fourier space, the extended 5N formalism (3.1) yields 


^ dN ^ ,, , 102AT f d^p d‘^N f d^p ^ 

Ck +^Ac,(k)+2^y +^j - p) 

at second order. 

The small scale perturbations are approximated to be of highly Gaussian distribu¬ 
tion, and thus, we have the following expressions, 

27r2 / H \ 2 

((5<y)k) = 0 , '^(k + k') , {d(t>]iS(j)]^'6(j)]^n) = 0 , (3.10) 

where we have defined k = |k|. As the large scale perturbation Ay is not stochastic, it has 
no correlation with 5(f) and itself. Using these relations and for ki ^ k, one obtains the 
asymmetric power spectrum 


(CkiCk,) 

27r2 (27r)^ 


-dN 
. dcj) 



2 


P^{k,v) 


(3.11) 















To apply Eq. (3.6), we can expand the power spectrum to the leading order of • (x — xq) 
as follows, 


Pf(A:,x) 


-ON 
- d(f) 


d^N 

dcgdcj) 


2 


Acs{t) 




2{d‘^ N / dcsdfj)) 


m I 
d<j> 


d^N 

dcs dip 


ACs{t) 


Ac<;(x,t) 


(3.12) 


As indicated in 
can rewrite the 
then derive 


As a result, the 

( 1 . 1 ). 


Eq. (3.8), Acs is only an overall shift of the background value. Therefore, we 
background sound speed parameter during inflation to be = cj + Ac* and 


ON d'^N _ dN 

d4> dcsdcp ^ dcp 


(3.13) 


asymmetric power spectrum can be further simplified as the same form of Eq. 


Moreover, the amplitude of the power asymmetry in this case is expressed as 


A{k) = ■ 

-(kLXis) 

d^N 

dcgdcf) 

rdN\-idcs . 

(3.14) 

Recall that, according to Eq. 

(2.15), 

the background dynamics of our 

model yields the 

following relation 

dN 

H 

1 

(3.15) 


d(j) 


Mpiy/2ecs 


Furthermore, by taking the partial derivative of the above formula with respect to Cg and 
considering the variation of the slow roll parameter e, we can get another useful relation 


^ =^_ (i + H) , 

dcsd(j) sJ ’ 


(3.16) 


where we have already applied the definition of slow roll parameters (2.11) and the following 
relation 


de e erj 

dCs Cs CgS ' 


(3.17) 


Finally, substituting (3.15) and (3.16) into (3.14), we obtain the amplitude of the power 
asymmetry of the following form 


A{k) = kixis (l + g ^ (l “ ^ • 


(3.18) 


Before moving to the next section, we would like to remark that in the EKC mechanism 
the asymmetry amplitude is related to the local non-Gaussianity via 


rlocal 
JNL 


5d^N //5iV\2 
6 / \ d(j) ) 


(3.19) 


Thus the observational constraints on primordial non-Gaussianity may result in the unnatural 
fine tuning of model parameters or a running behavior of if the power asymmetry is due 
to the inflaton field[36, 41]. In our approach the asymmetry amplitude (3.14) is not related 
to the primordial local non-Gaussianity, and thus, there is no necessity to require any fine 
tuning or a running 
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4 Specific example and confronting observations 


In this section, we use specific implementation of the multi-speed model to present how our 
mechanism provides a successful explanation for the hemispherical power asymmetry and use 
the observational data to constraint the model. 

Notice that the constraint of Planck on the primordial non-Gaussianity of the DBI model 
only gives = 2.6 ± 61.6 and Cg > 0.069 [60], thus there is a large parameter range that 

we could explore in our case. For example, if we demand ~ 1 (i.e. /ati, <C 1 ) and r]/s <C 1, 
Eq. (3.18) reduces to 

|A(m = (4.1) 

with = kL'r’isB{t) sin6. This result has been fully discussed in the previous work [29]. 

Now let us take a more specific example based on Ref. [54] which indicates the case 
with r/ ~ s and <C 1. As we have pointed out, there are two undefined functions in the 
original Lagrangian, the warping factor f{4>,x) the potential V{(j),x)- If we consider an 
AdS warping throat, the warping factor can be expressed as 

/(</>,X) =/(0)+/(x) = ^, (4.2) 


while the IR type potential takes the following form 

=Vo- _ (4 3) 

Under the relativistic limit of inflaton with small sound speed and the requirement \Ht\ S> 1 
during inflation, an approximate solution of Eq. (2.9) is given by [54, 55] 


(j) = 


y^(i _ 
t V 2mH'^ 



(4.4) 


where the beginning of inflation is set at t —)• — oo and higher order terms are suppressed by 
By Eq. (2.3) the leading order of the sound speed parameter takes the form of 






2>H 

m?t 


(4.5) 


by 


Afterwards, one can calculate the slow roll parameters at leading order, which are given 


Xm? 


e = 






_ e 3 I ^ 

^ ^ Hcl ^ 




) ■ 


(4.6) 


This indicates that tj/s 
(3.18) yields 


3 in the relativistic limit. Thus the asymmetry amplitude in Eq. 


|A(fc)|^^-7^R(t)sinR 


x(4) Cs(4)^ 


(4.7) 











From Eq. (4.7), we can see that there exists a fc-dependent behavior of the asymmetry 
amplitude which originates from the evolution of the x held via the sound speed during 
inflation. The amplitude of the power spectrum is fixed when the mode k exits the Hubble 
radius, i.e. k < aH. The perturbation modes of different scales correspond to different 
Hubble-exit time while the evolution of x well as Cg gives different field values at 
Thus this can be used to provide a natural explanation of the observational fact that A{k) is 
significant at cosmological scale but turns to be negligible at the Mpc scale. 

For the amplitude of the super-horizon perturbation Ay, it is unclear to us whether 
this amplitude is enhanced or not. Since there is no observational evidence to imply any 
enhancement of perturbation modes at super-Hubble scales, it is reasonable to assume that 
Ax has the same amplitude with 5x and 5(f). Hence, we can approximately take kLXisB{t) ~ 

^ = y/2eCsP^^'^Mpi ~ AIpi- Then Eq. (4.7) becomes 


\A{k)\ ~ 


l-Ug Pl^'^pMpi sin 9 
2 x{tk)cg{tkf/'^ 


(4.8) 


According to the latest Planck data [60, 61], we learn that Ug = 0.968 ±0.006, Cg > 0.069 and 
ln(10^*^P^) = 3.062 ± 0.029 at the 68% confidence level. In the following, we set Cg{tg) = 0.1 
at the Planck pivot scale kg = 0.05Mpc“^. Then by applying Eq. (2.4) we immediately get 
= —32, which is consistent with the constraint of the Planck data, which is = 

2.6 ± 61.6 [60]. Making use of Eq. (4.5) and horizon exit condition a(tk)H{tk) = k, we get 
the /c-dependent sound speed, which is expressed as 


Cg{tk) 


Cgitgf 



m^Cgftg) 

-— In 

3772 



(4.9) 


Further, by taking psin0 ~ 1, we obtain 


\A{k)\ ~ 


18.5 X 10-^Mpi 
x{tk) 



rn^Cgftg) 



(4.10) 


For the large scale mode with7 < 60, the CMB observation yields A = 0.072 ± 0.022. 
Considering k ~ l/xig ~ 7.5 x 10“®7, we take A = 0.072 at ka ~ 0.0045Mpc“^, which exits 
the Hubble radius at Therefore, we obtain 


X(4,) - 4.38 X 10 ^ pMpisin6» . (4.11) 

However, for the small scale mode with ki, ~ lMpc~^, which exits the Hubble radius at tk,^, 
the quasar observation indicates A ~ 0.0153. Correspondingly, we derive 


X(4j - 4.25 X 10 ^ pMpi sine . 


(4.12) 


Since the large scale mode exits the Hubble radius much earlier than the small scale mode, 
the calculation above suggests that during inflation the x held evolves from large values to 
small ones, which happens to be the picture of large field inflation. 

To realize this picture, we turn to study the evolution of x- Since it is a light field during 
inflation, we can neglect the effect of potential term, i.e. V (x) — 0. Thus the equation of 
motion of this canonical field is x ± 377x = 0; which yields 

X{t) = C - (4.13) 
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for a nearly constant Hubble constant. Here C and D are the integral constants and can be 
determined by Eq.(4.11) and Eq.(4.12). With this solution and horizon exit condition, we 
have xi^k) = C — D'k~^. Then the asymmetry amplitude in Eq.(4.10) is determined, whose 
scale-dependence is shown in the left panel of Fig.l. Furthermore, we can also demonstrate 
the variation in CMB temperature power spectrum by = f where Axe(k) 

is transfer function. Then 

ACr !2A{k)P^AUk)dk/k 

f P(^A'^^{k)dk/k 

Considering that the transfer function is sharply peaked at i = kxis, we plot the variation of 
the TT spectrum in the right panel of Fig.l. 




Figure 1. Left: The asymmetry amplitude A varies with the comoving wave number k. Right: The 
fractional changes ACj"^ jCj"^ in the CMB TT spectrum. The shadow and red error bar correspond 
to the observational constraints at different scales, which fix the undetermined coefficients. 


5 Conclusion 

The hemispherical asymmetry of the CMB temperature fluctuations has drawn a lot of at¬ 
tention of cosmologists. This anomaly, if confirmed, will be an important result which is 
believed to connect with some unknown new physics. Thus it is interesting to explore plau¬ 
sible theoretical explanations. In Ref. [29], it was found that a statistically inhomogeneous 
sound speed parameter can lead to the spatial variation of the primordial spectrum and a 
rough calculation of the asymmetry amplitude was provided. 

In the present paper, we firstly review the model of multi-speed inflation and the basic 
ideas of 5N formalism. Then, we generalize the 6N formalism to include consistently the 
contribution of the spatially varying sound speed. The modulation of Cg can be caused by a 
perturbation mode of very long wavelength seeded by a light field which is not necessary to 
contribute to the background evolution as well as the generation of curvature perturbations 
during inflation. And then in the framework of the multi-speed inflation, we give a system¬ 
atic calculation of the asymmetry amplitude. We find that our approach is free from the 
constraints of the local non-Gaussianity. For some particular choices of model parameters, we 
can deduce back to the result derived in [29]. However, for a more generic case, the behavior 
of the spatial variation presents significantly differences from the previous result. We point 
out that our solution may provide a plausible explanation of the scale dependence of the 
asymmetry amplitude, which accommodates the observations of both CMB and quasar. 
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Finally we would like to remark that the generalized 6N formalism, which has taken 
into account the modulation of primordial sound speed parameter, can be applied to any 
early universe models involving nonstandard kinetic terms. Therefore, besides the specific 
multi-speed model we considered in the present paper, it is also interesting to further study 
inflation models with other k-essence fields. 
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